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PREFACE 


Panel method code ANTARES is described in this report. The code computes sub- 
sonic wall interference effects in a wind tunnel with rectangular cross-section. Different 
types of wall boundary conditions may be specified by the user. Blockage effects of a wind 
tunnel model are represented by point doublets. Lifting effects of a wind tunnel model are 
represented by line doublets. Compressibility effects are modeled by applying the Prandtl- 
Glauert transformation. Wall interference correction calculations were successfully verified 
by comparing panel method code solutions with corresponding classical solutions. A mod- 
ified version of panel method code ANTARES is used to compute perturbation velocity 
database files for the real-time wall interference correction system of the NASA Ames lift 
Transonic Wind Tunnel. 

I hope that the promising results obtained during the development of panel method 
code ANTARES will benefit attempts to improve data quality and efficiency of wind tunnel 
tests in the NASA Ames lift Transonic Wind Tunnel. 

I want to thank Alan Boone, Mike George, and Don Nickison of NASA and Doug Pena 
of SVERDRUP Technology for their critical and constructive review of the manuscript. 
Their careful and competent proofreading and checking has greatly improved the quality 
of this report. 
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May 2000 
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ABSTRACT 


Panel method code ANTARES was developed to compute wall interference corrections 
in a rectangular wind tunnel. The code uses point doublets to represent blockage effects and 
line doublets to represent lifting effects of a wind tunnel model. Subsonic compressibility 
effects are modeled by applying the Prandtl-Glauert transformation. The closed wall, open 
jet, or perforated wall boundary condition may be assigned to a wall panel centroid. The 
tunnel walls can be represented by using up to 8000 panels. The accuracy of panel method 
code ANTARES was successfully investigated by comparing solutions for the closed wall 
and open jet boundary condition with corresponding Method of Images solutions. Fourier 
transform solutions of a two-dimensional wind tunnel flow field were used to check the 
application of the perforated wall boundary condition. Studies showed that the accuracy 
of panel method code ANTARES can be improved by increasing the total number of wall 
panels in the circumferential direction. It was also shown that the accuracy decreases with 
increasing free-stream Mach number of the wind tunnel flow field. 
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CHAPTER 1 


INTRODUCTION 

The calculation of wind tunnel wall interference corrections in a perforated or slotted 
wall wind tunnel with rectangular cross-section has always been a challenge due to the 
complexity of the flow field in the vicinity of the test section walls. Numerical solutions 
of the wall interference flow field have to be computed by using a sufficiently accurate 
description of the wall boundary conditions. In addition, compressibility effects have to 
be included in the wall interference calculation as perforated or slotted wall wind tunnels 
are frequently operated at high subsonic Mach numbers. 

Presently, a real-time wall interference correction system is being developed for the 
slotted wall test section of the NASA Ames lift Transonic Wind Tunnel (TWT). This 
system is similar to a system that was successfully implemented in the NASA Ames 12ft 
Pressure Wind Tunnel [1]. The real-time correction system uses a singularity represen- 
tation of the test article. The system also requires precomputed solutions of the wall 
interference flow field in the test section and of the pressure coefficient on selected wall 
pressure ports. Therefore, it became necessary to develop a three-dimensional flow field 
solver, i.e. panel method code ANTARES, that allows for the calculation of the wind 
tunnel and wall interference flow field of a singularity in a perforated or slotted wall wind 
tunnel with rectangular cross-section. 

Different numerical techniques are available to compute wall interference in a perfo- 
rated or slotted wall wind tunnel. After careful review of the existing literature [2] the 
author decided to use Keller’s panel method algorithm [3] as a basis for the development 
of panel method code ANTARES. 

Keller [3] developed an incompressible panel method algorithm that computes the 
upwash velocity (angle of attack correction) caused by line doublets in a wind tunnel 
with rectangular cross-section. Keller’s description of the panel method algorithm is 
exceptionally clear and complete. His algorithm allows the user to assign one of six different 
types of boundary conditions to each wall panel. The following boundary conditions may be 
specified: (1) closed wall; (2) open jet; (3) perforated wall; (4) ideal slotted wall (integrated 
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form); (5) ideal slotted wall (differentiated form); (6) slotted wall including viscosity in 
slots. Slopes of the source distribution strength of each wall panel are unknowns in the 
linear system of Keller’s algorithm. The upwash velocity caused by the effect of the test 
section walls is computed as soon as the solution of this linear system, i.e. the solution of 
the boundary value problem, is obtained. 

Panel method code ANTARES uses a modified version of Keller’s algorithm to set up 
the linear system that describes the wall boundary condition at each wall panel centroid. 
Several additions and extensions were made to Keller’s algorithm so that panel method 
code ANTARES may be used to compute blockage and angle of attack corrections in a 
subsonic wind tunnel flow field. Point doublets were chosen for the calculation of blockage 
corrections and the Prandtl - Glauert transformation was selected to model compressibility 
effects. It was also decided to include features in panel method code ANTARES that allow 
for the calculation of semispan model wall interference effects. 

This report describes basic assumptions and elements of panel method code ANTARES. 
At first, Keller’s algorithm [3] is discussed. Then, the modeling of compressibility effects is 
explained in great detail. A description of the structure of panel method code ANTARES 
follows. Finally, wall interference calculation results for different types of wall boundary 
conditions are compared with available classical solutions. 
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CHAPTER 2 


PANEL METHOD CODE ANTARES 
2.1 Wall Boundary Condition Description 


Panel method code ANTARES was developed to compute the wind tunnel and wall 
interference flow field in a wind tunnel with rectangular cross-section. ANTARES uses 
Keller's algorithm [3] to set up the matrix equation that allows the user to describe six 
different types of wall boundary conditions. 

In incompressible flow, the wall boundary 7 conditions may be expressed by using the 
perturbation velocity potential 0 of the wind tunnel flow field [3]. Then, we get : 


, d 0 

C\ • 0 + C2 • — ^ 

a x 


+ 


C3 • 


d <i> 

d n 


+ 


d 2 

° 4 d x d n 


0 ( 1 ) 


where .x is the streamwise coordinate, n is the normal vector on the wind tunnel wall, and 
ci,C 2 ,C 3 ,C 4 are coefficients that describe the different types of wall boundary conditions. 
The following table lists values of these coefficients for six different types of boundary 
conditions: 


Table 1 : Wall Boundary Condition Coefficients (from Ref. [3]) 


Type of Boundary Condition 

Cl 

C2 

C3 

c 4 

Closed Wall 

0 

0 

l 

0 

Open Jet 

0 

i 

0 

0 

Perforated Wall 

0 

l 

l/R 

0 

Ideal Slotted Wall (Integrated Form) 

1 

0 

l 

0 

Ideal Slotted Wall (Differentiated Form) 

0 

1 

dl/dx 

/ 

Slotted Wall Including Viscosity in Slots 

0 

1 

dl/dx + l/R 
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The total perturbation velocity potential <f> may be written as the sum of the potential 
(f> m caused by the test article in free-air and the potential <p w caused by the wind tunnel 


3 
























walls, i.e. the wail interference potential. We get : 
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( 2 ) 


Combining Eqs. (1) and (2), we get : 
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In general, a panel method code is constructed such that boundary conditions are only 
fulfilled at each wall panel centroid. Each panel may be considered as an infinitesimal wall 
element with its own boundary condition characteristics. Thus, wall boundary condition 
coefficients cj, • • • , c 4 axe a function of each wall panel centroid “i” . Using Eq. (3), the 
boundary c edition at each panel centroid may be written as : 
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where [4> m ]i is the model potential and [<f> w ]i is the wall potential at the panel centroid 
“i”. The model potential 4> m of a point doublet and corresponding derivatives are given in 
App. 1. The model potential <f> m of a line doublet and corresponding derivatives are given 
in App. 2. It remains to determine the wall potential 4> w . 

Keller [3] selected a wall panel element in his algorithm that may be used to express 
the wall potential [<t> w \i at each panel centroid as a function of the source strength slope 
fij of all “IV” panel elements describing the wind tunnel. Then, we get [3]: 
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3 = 1 

where [<j >*\ ,, is the perturbation potential per unit source strength slope at panel centroid 
“i” due to panel “j”, and nj is the source strength slope of panel “y”. The perturbation 
potential <t>* and its derivatives are given ; .d Ref. [3]. Combining Eqs. (4) and (5), we get : 
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Equation (6) has to be applied to each wall panel centroid “i” (i.e. i = 

Thus, a linear system for the unknown source strength slope fij of each wall panel "j" (i.e. 
/i!, • • • ,^n) is obtained. This linear system may be written as : 
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hi = - C X (0 ’ - <*(*)* ~ C 3(0 

Finally, after the solution of the source strength slope of each wall panel is ob- 
tained by applying a linear system solver to Eq. (7a), it is possible to compute the axial 
perturbation velocity at a wall panel centroid “i” in the wind tunnel flow field. Combining 
Eqs. (2) and (5) and differentiating the result in the streamwise direction we get : 
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The pressure coefficient at a wall panel centroid may be computed as : 
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Similarly, perturbation velocity components of the wall interference flow field at a flow 
field point “fc” may be written as : 
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Using Eq. (9a), we get the following equation for the blockage factor at a flow field 
point u k n in the test section : 


e(k) = 




N r 


Ur 


t-'-E 

u no _ 


d <f>* 


j=\ «■ 




(10a) 




Similarly, assuming that the lift force of the test article points in the positive z- 
direction (see Fig. 1) and using Eq. (9c), we get for the angle of attack correction in 
degrees : 
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2.2 Compressibility Effects 

The calculation of wind tunnel and wall interference effects of a point source, point 
doublet, or line doublet located in a rectangular wind tunnel requires the application 
of some sort of compressibility correction if tests sure conducted at a high subsonic Mach 
number. Unfortunately, Keller’s [3] panel method code algorithm is limited to incompress- 
ible flow. Therefore, it became necessary to include the Prandtl-Glauert compressibility 
transformation in panel method code ANTARES. 

In general, the Prandtl-Glauert transformation has to be applied to solve the subsonic 
potential equation in a subsonic flow field. A detailed description of the Prandtl-Glauert 
transformation may be found in Ref. [4] . The subsonic potential equation is given as : 


fi . 52 * + 
P dx 2 

d 2 cj> 
dy 2 + 
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II 

CO <0 

(lla) 


= 1 - 

M 2 

1± oo 

(114) 


where <f> is the perturbation velocity potential, x,y,z are cartesian coordinates of the flow 
field, and M is the free-stream Mach number of the subsonic flow field. 


6 



The subsonic potential equation is a partial differential equation that has to be solved 
for the given wind tunnel wall boundary conditions and for the selected singularity (i.e. 
point source, point doublet, or line doublet) that represents the test article for the purpose 
of wail interference calculations. This solution is basically obtained in three steps: 

(1) The subsonic potential equation and corresponding boundary conditions axe trans- 
formed by applying the Prandtl— Glauert transformation. 

(2) The transformed subsonic potential equation is a Laplace equation; the trans- 
formed boundary value problem is solved numerically by using a panel method code. 

(3) The solution of the original subsonic potential equation is obtained by reversing 
the Prandtl-Glauert transformation; perturbation velocity components computed in the 
transformed coordinate system are transformed back to the original coordinate system. 
Variable Transformation 

The Prandtl-Glauert transformation maps variables from an original coordinate sys- 
tem x,y,z to a transformed coordinate system x' ,y',z' . The transformation of the coordi- 
nates x,y,z , of the free-stream velocity Uoc, and of the perturbation potential <j> is given 
by the following equations (see Eq. (4-10) in Ref. [4]) : 
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~ X 
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In the case of a subsonic wind tunnel flow field, it is also necessary to apply the 
Prandtl-Glauert transformation to variables that are used to describe wind tunnel wall 
boundary conditions. Thus, it will be shown below that the restriction parameter R 
(perforated wall boundary condition) has to be transformed using the following equation : 


t 



R = -f- (12/) 


Similarly, it will be shown below that the slot parameter l (slotted wall boundary 
condition) has to be transformed using the following equation : 


v = i -(3 ( 12 $) 


The normal vector n at the wind tunnel wall is located in the y-z-plane. Thus, we get: 
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The transformed normal vector n' may be written as : 
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Combining Eqs. (12b), (12c), with Eq. (12i), we get : 
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Comparing Eq. (12j) with Eq. (12h), we get : 
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Wall interference calculations may use singularities to represent the blockage and 
lifting effect of a test article. The Prandtl-Glauert transformation has to be applied to the 
location of the singularities. The strength a of the singularities is connected with the test 
article geometry and has to be transformed as well. The transformed singularity strength 
<r' may be expressed as : 
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where F(/3 ) is a transformation function. The value of this transformation function is 
rigorously derived in App. 3 for different types of singularities (see also Table 7 on p.54). 
Subsonic Potential Equation Transformation 

The subsonic potential equation given in Eq. (11a) may now be transformed. Differ- 
entiating Eqs. (12a), (12b), and (12c), we get : 
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Applying the chain rule to Eq. (12e) twice and using Eqs. (13a), (13b), (13c), we get : 
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Finally, the transformed subsonic potential equation is obtained by combining Eqs. (11a), 


(14a), (14b), and (14c) . We get : 
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In the next step it is necessary to apply the Prandtl-Glauert transformation to the 
wall boundary conditions. 
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Wall Boundary Condition Transformation 

Keller [3] discusses six different types of wall boundary conditions. The Prandtl- 
Glauert transformation has to be applied to each type before a solution of the subsonic 
flow field can be found. 

Type 1 / Closed Wall : The closed wall boundary condition may be written as : 
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From Eq. (12h) we know that N x = 0 on the wind tunnel wall. Thus, we get for Eq. (16b) : 
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Similarly, the transformed normal derivative operator may be written as : 
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Combining Eqs. (13b), (13c), with Eq. (16e), we get : 
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From Eq. (12k), we know that N' y = N y and N' z = N z . Thus, using Eq. (16c), we get for 
Eq. (16e) : 
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Combining Eq. (16a) with Eqs. (12e) and (16g), we get : 
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Type 2 / Open Jet : The open jet boundary condition may be written as : 
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Applying the chain rule to Eq. (17a), we get : 
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Combining Eqs. (12e), (13a) with Eq. (17b), we get : 
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Type 3 / Perforated Wall : The perforated wall boundary condition may be written 
as : 
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Applying the chain rule to Eq. (18a) and using Eq. (16g), we get : 
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Combining Eqs. (12e), (13a), with Eq. (18b), we get : 
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It is possible to define a transformed restriction parameter R! = R/ 3 ■ Then, the 
transformed boundary condition for a perforated wall becomes : 


d <}>' 1 d <{>' 

d 7' + R' dri 


0 (18e) 


Type 4 / Ideal Slotted Wall (integrated form) : The ideal slotted wall boundary 
condition in integrated form may be written as : 


<t> + l- = 0 (19a) 

d n 


Using Eq. (16g) in Eq. (19a), we get : 


4> + l ■ 3 


d <j> 
d n' 


= 0 


(196) 


Combining Eq. (12e) with Eq. (19b), we get : 




3 2 


d n' 


3 2 


4>' + l • 3 • 


d <f>' 
d n' 


= 0 (19c) 


<t>' + { l ’ 0 ) 


d <y 

d n' 


= 0 


(19cQ 
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A transformed slot parameter V = l • 0 is introduced in Eq. (19d). Then, comparing 
Eqs. (19a) and (19d), the transformed boundary condition for an ideal slotted wall in 
integrated form becomes : 


=> 0' + /'• 


d 6' 

d n' 


0 (19e) 


Type 5 / Ideal Slotted Wall (differentiated form) : The ideal slotted wall boundary 
condition in differentiated form may be written as : 


d <i> d l d 0 i d 2 <t> 
d x d x d n + d x d n 


0 (20a) 


Applying the chain rule to Eq. (20a) and using Eq. (16g), we get : 

d x' d 0 d x' d l d 0 dx' d 2 0 

— 4. g . _ 4. l 3 — 

d x d x' d x d x 1 d n' d x d x' d n' 


0 (206) 


Combining Eqs. (12e), (13a), with Eq. (20b), we get : 


<,(*'//>’) + 4±. e . 8i*J0‘) +l . B . £j*Li£) . o (20c) 


d x 1 


d x 7 


d n' 


a * + 2LL£I.44 + (i ., ) 


d x' d n' 

d 2 <j>' 


= 0 


(20 d) 


dx' dx' dri v " y dx' dri 

A transformed slot parameter /' = /•/ 3 is introduced in Eq. (20d). Then, the trans- 
formed boundary condition for an ideal slotted wall in differentiated form becomes : 


d£ a_e_ _a_f_ a* y 

d x' d x' d n' d x' d n' 


0 (20e) 


Type 6 / Slotted Wall Including Viscosity : The slotted wall boundary condition 
including viscosity in slots may be written as : 
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Applying the chain rule to Eq. (21a) and using Eq. (16g), we get : 

d x' d <{> t d x^ d l j 1_ R d ^ f d x' ^ d 2 (/> 

3 x d x' + d x d x 1 R d n' d x d x' d n' 


Combining Eqs. (12e), (13a), with Eq. (21b), we get : 


d x‘ 


d I , 1 1 a a ( 6' / p ) , f .a 2 ( <t>' / /? 2 ) 
~dV R J' dn' + P ' dx> dn' 


0 ( 216 ) 


0 (21c) 


d(l-0) 


- + (/■/?). 


d 2 

d x' d to! 


= 0 (21d) 


d x' [ d x' (R/ 0) \ dn' K ^ J d x' d n' 

Now a transformed restriction parameter R' = R/fi and a transformed slot parameter 
V = l-P are introduced. Then, comparing Eqs. (21a) and (21d), the transformed boundary 
condition for a slotted wall including viscosity in slots becomes : 



Reversal of the Prandtl-Glauert Transformation 

The solution of the transformed subsonic potential equation (Eq. (15)) may be com- 
puted by using the panel method code algorithm developed by Keller [3]. After the solution 
is found, it is only necessary to reverse the Prandtl-Glauert transformation to obtain the 
solution of the subsonic wind tunnel flow field. Perturbation velocities computed in the 
transformed coordinate system have to be transformed back to the original coordinate 
system. Thus, applying the chain rule to Eq. (12e), we get for the perturbation velocity 
components in the subsonic flow field : 

d 4> 1 d x' d <f>' 

U = d x = ~fP JFx dUF 


14 


(22a) 




V = 

d d> 
d y 

= 

1 

P 2 

dy' 

d y 

d4>' 
d y' 


(226) 

W — 

d <t> 

d z 

= 

1 

Q 2 

d z' 
d z 

dp 

d z' 


(22c) 

Finally, combining Eqs. (13a), (13b), (13c) with Eqs. (22a), (22b), (22c), 
perturbation velocity components in the subsonic flow field : 

we get the 


U — 

i 

,d 2 

d <t>' 

d x* 

— 

u f 

P 2 

(23a) 




V — 

i 

T' 

d <t>' 
d y' 

— 

v' 

(3 

(236) 




w = 

i 

T‘ 

d 4>' 
d z' 

= 

w l 

0 

(23c) 




where u', i/, w' are the perturbation velocity components in the transformed coordinate 
system x' , y', z' . 

In general, the application of the Prandtl-Glauert transformation in a panel method 
code may be summarized as follows : 

(1) In a first step, the original coordinate system (r,y,z) maps to a transformed 
coordinate system (x 1 ,y', z see Eqs. (12a), (12b), (12c)). Restriction parameter R and slot 
parameter l have to be transformed by using Eqs. ( 1 2f ) , (12g). Strengths a of singularities 
representing the test article have to be transformed by using Eq. (121) and Table 7 on 
p.54 . 

(2) In the next step, the solution of the transformed subsonic potential equation is 
found by using the algorithm given by Keller [3]. 

(3) Finally, perturbation velocities in the subsonic wind tunnel flow field are computed 
by reversing the Prandtl-Glauert transformation (Eqs. (23a), (23b), (23c)). The following 
flow chart summarizes different solution steps : 
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Subsonic Wind Tunnel and Wall Interference Flow Field 

II 

Subsonic Potential Equation (Eqs. ( 11 a), (lib)) 

il 

Prandtl-Glauert Transformation (Eqs. ( 12 a), ( 12 b), ( 12 c), (12f), (12g), (121)) 

Laplace Equation (Eq. (15)) 

II 

Panel Method Code Solution ( Keller’s Algorithm, Ref. [3]) 

II 

Reversal of Prandtl-Glauert Transformation (Eqs. (23a), (23b), (23c)) 

II 

Solution of Wind Tunnel and Wall Interference Flow Field 


2.3 Panel Method Code Structure 


Panel method code AN TARES computes the wall interference flow field in a wind 
tunnel with rectangular cross-section. ANTARES allows the user to represent the test 
article by using two types of singularities. Point doublets represent blockage effects of the 
test article. Line doublets represent lifting effects of the test article. Currently, ANTARES 
allows the user to select three different types of wall boundary conditions (closed wall, open 
jet, or perforated wall). Other boundary conditions listed in Table 1 on p.3 may easily be 
implemented as the code computes all derivatives required for the application of Keller’s 
algorithm [3]. The internal structure of panel method code ANTARES can be divided into 
three parts (see Fig. 2 ) : 

(i) Preparation of Linear System: In the first part, the linear system describing 
the wall boundary conditions is prepared. At first, the wall boundary condition type of 
each wall panel centroid “i” has to be specified by selecting coefficients c\ (i) to c 4 (z) . These 
coefficients are required for the calculation of matrix coefficients a t} and the right hand side 
vector bi of the linear system (see Eqs. (7b), (7c)). The Prandtl-Glauert transformation 
has to be applied to the restriction parameter “/?” and to the slot parameter if they 
are used to compute coefficients 03 ( 1 ) and c 4 (i). The reciprocal of restriction parameter 
“JF is used to compute coefficient 03 ( 1 ). Therefore, panel method code ANTARES checks 
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if the absolute value of “R” is greater than 0.0001 . 

In the next step the panel model of the wind tunnel walls has to be prepared. Equal 
panel spacing is selected in the circumferential direction of the test section. Cosine panel 
spacing is selected in the streamwise direction of the test section with smaller panels 
near the location of the test article. Experience has shown that a total number of 60 
panels in the streamwise direction is sufficient as Keller’s algorithm uses a wall panel type 
that linearly varies the source strength over each panel in the steamwise direction [3]. 
However, a relatively large number of panels is required in the circumferential direction 
of the test section as Keller’s panel type keeps the source strength constant over each 
panel in the circumferential direction. Thus, a total number of 4800 panels (80 panels 
in the circumferential direction and 60 panels in the streamwise direction) is selected to 
represent the test section for a fullspan model configuration (see Fig. 4). A total number 
of 7200 panels (120 panels in the circumferential direction and 60 panels in the streamwise 
direction) is selected to represent the test section for a semispan model configuration 
(see Fig. 8). For a fullspan model configuration, the length of the paneled test section 
geometry is about eight tunnel widths. Finally, after the Prandtl-Glauert transformation 
is applied to the wall panel coordinates, it is possible to compute matrix coefficients a,j 
using Eq. (7b). 

Type and location of singularities representing the test article have to be specified 
as well so that the right hand side vector bi of the linear system can be computed. The 
Prandtl-Glauert transformation has to be applied to singularity coordinates and the sin- 
gularity strength. Afterwards it is possible to apply Eq. (7c) and compute the right hand 
side vector b{ of the linear system given in Eq. (7a). 

(ii) Solution of Linear System: In this part a linear system solver is used to 
solve the boundary value problem. Numerical studies have shown that a standard off- 
the-shelve LU-Decomposition algorithm [5] may be used to solve the linear system. It 
is recommended to use “DOUBLE PRECISION” arithmetic in an effort to avoid finite 
digit arithmetic problems. ANTARES stores all matrix coefficients in arrays and therefore 
storage requirements can be large. For example, the storage of coefficients aij for 4800 
panels requires approximately 184 [Mbytes]. ANTARES takes advantage of symmetry in 
the case of the test section panel model for a semispan model configuration. This reduces 
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the number of unknowns in the linear system from 7200 to 3600. 

(iii) Calculation of Wall Interference Flow Field: After the solution of the linear 
system defined by Eq. (7a) is found wall interference perturbation velocity components 
may easily be computed by adding all velocity contributions of wall panels using Eqs. (9a), 
(9b), (9c). Finally, the Prandtl-Glauert transformation has to be reversed using Eqs. (23a), 
(23b), (23c) so that perturbation velocity components in the compressible subsonic wall 
interference flow field may be obtained. 

The real-time wall interference correction system currently being developed for the 
NASA Ames lift TWT requires the calculation of the pressure coefficient at each wall 
panel centroid. ANTARES uses Eq. (8b) to compute this pressure coefficient. 

2.4 Closed Wall Boundary Condition 

The accuracy of panel method code ANTARES is investigated in this chapter. There- 
fore, solutions of the wall interference calculation for a closed wall rectangular wind tunnel 
are compared with corresponding Method of Images solutions provided in App. 4 and 
App. 5. In addition, solutions of the wind tunnel flow field on the wind tunnel walls are 
investigated. 

In general, panel method code ANTARES may be used to compute the wall inter- 
ference flow field caused by a point doublet or line doublet in a closed wall wind tunnel. 
The user only has to specify the panel representation of the test section geometry and the 
location of the selected point doublet or line doublet inside of the test section. Coeffi- 
cients related to the wall boundary condition type on each of the four test section walls 
(see nomenclature in Fig. 1) have to be specified as well. These coefficients are defined in 
Eq. (1). The table below lists the coefficients for the closed wall boundary condition : 


Table 2 : Wall Boundary Condition Coefficients (Closed Wall) 


Wall No. 

Panel Index Range 

ci(0 

c 2 {i) 

C 3 (0 

C*(0 

1 

1 < i < 1200 

0 

0 

*1 

i 

0 

2 



0 

1 

0 

3 

2401 < i < 3600 

0 


1 

0 

4 

3601 < i < 4800 

0 

0 

1 
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Two different wind tunnel cross-sections are selected for the study of the accuracy of 
panel method code ANTARES. Compressibility effects are investigated as well by comput- 
ing wall interference for two different Mach numbers (0.00 and 0.80). 

Wind Tunnel Cross-Section No.l : The first wind tunnel used for the accuracy 
study has a square cross-section (see Fig. 3). The panel representation of the test sec- 
tion boundary is depicted in Fig. 4. The wall panel index range of each wall is given 
in Table 2 . A total number of 4800 panels was selected for the panel model (60 panels 
in the streamwise direction and 80 panels in the circumferential direction). Cosine panel 
spacing was selected in the streamwise direction and equal panel spacing in the circum- 
ferential direction. The total length of the paneled test section was chosen to be eight 
tunnel widths. The calculation of wall interference effects caused by a fullspan model is 
simulated. Therefore, a single point doublet or line doublet is placed inside of the t unn el 
to represent blockage or lifting effects of the wind tunnel model (see Fig. 3). The selected 
tunnel width and height is 10 [ft]. For the present study it was decided to compute the 
wall interference flow field on the tunnel centerline (see Line 1 in Fig. 3). The wind tunnel 
flow field on selected locations on the wind tunnel wall is computed as well (see Row 1 and 
Row 2 in Fig. 3). 

Figures 5a, 5b, and 5c depict the result of the wall interference calculation for the 
selected unit point doublet representing test article blockage. The dimensionless wall 
interference perturbation velocity components are depicted as a function of the t unn el 
centerline coordinate. The overall agreement between panel method code solution (abbre- 
viated P.M.) and Method of Images solution (abbreviated M.I.) for both Mach n um bers 
is excellent. It appears that the difference between the panel method code solution and 
the Method of Images solution (exact solution) is small for the y and z component of the 
wall interference velocity vector (see Fig. 5b and 5c). Differences between panel method 
code solution and the Method of Images solution (exact solution) increase with increasing 
Mach number only for the axial interference velocity component (see Fig. 5a). 

Figures 5d and 5e compare the wind tunnel flow field solution on Row 1 and Row 2 for 
the selected unit point doublet. Again, the overall agreement between the panel method 
code solution and the Method of Images solution for the two selected Mach numbers is 
good. 
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Figures 6a, 6b, and 6c depict the result of the wall interference calculation for a unit 
line doublet representing test article lifting effects. The line doublet orientation angle is se- 
lected to be 0.0°. Again, dimensionless wall interference perturbation velocity components 
are shown. The agreement between panel method code solution and Method of Images 
solution for all three components of the wall interference velocity vector and for both Mach 
numbers is good. Differences between the panel method code solution and the Method of 
Images solution (exact solution) are smaller than in the case of the unit point doublet. 

Figures 6d and 6e compare the wind tunnel flow field solution on Row 1 and Row 2 
for the selected unit line doublet. Again, the agreement between the panel method code 
solution and the Method of Images solution for the two selected Mach numbers is good. 

Wind Tunnel Cross-Section No. 2 : The second wind tunnel geometry used for 
the present accuracy study has a rectangular cross-section (see Fig. 7). The corresponding 
panel representation of the test section is depicted in Fig. 8. A total number of 7200 panels 
was selected for the panel model (60 panels in the streamwise direction and 120 panels in the 
circ umf erential direction). Cosine panel spacing was selected in the streamwise direction 
and equal panel spacing in the circumferential direction. The total length of the paneled 
test section was chosen to be eight tunnel widths. The calculation of wall interference 
effects of a semispan model configuration is simulated in this example assuming that the 
reflection plane is located halfway between Wall 3 and Wall 4. A point doublet or line 
doublet and its corresponding mirror image are placed inside of the tunnel to represent 
lifting or blockage effects of a semispan wind tunnel model (see Fig. 7). The selected tunnel 
width is 10 [ft]. Wall interference of a semispan model mounted on the floor of a 10 [ft] 
by 10 [ft] test section is being computed (see Line 2 in Fig. 7). Therefore, the tunnel 
height has to be chosen 20 [ft]. The wind tunnel flow field on selected locations on the 
wind tunnel wall is computed as well (see Row 1 and Row 2 in Fig. 7). 

Figures 9a, 9b, and 9c depict the result of the wall interference calculation for the 
selected two unit point doublets representing test article blockage. Dimensionless wall in- 
terference perturbation velocity components are depicted. The overall agreement between 
panel method code solution and Method of Images solution for both Mach numbers is 
excellent. 

Figures 9d and 9e compare the wind tunnel flow field solution on Row 1 and Row 2 
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for the selected two unit point doublets. Again, the overall agreement between the panel 
method code solution and the Method of Images solution for the two selected Mach n um bers 
is good. 

Figures 10a, 10b, and 10c depict the result of the wall interference calculation for the 
selected two unit line doublets representing lifting effects of the semispan model. The line 
doublet orientation angle is selected to be 90.0° and therefore the lift force of the semispan 
model points in the positive y-coordinate direction (see Fig. 7). The dimensionless wall 
interference perturbation velocity components are depicted. The overall agreement be- 
tween panel method code solution and Method of Images solution for both Mach n um bers 
is good. 

Figures lOd and lOe compare the wind tunnel flow field solution on Row 1 and Row 2 
for the selected two unit line doublets. Again, the overall agreement is good. 

2.5 Open Jet Boundary Condition 

Panel method code ANTARES allows the user to compute the wind tunnel wall in- 
terference flow field of a test article in an open jet. Method of Images solutions of a point 
doublet or line doublet in a rectangular open jet are used to study the accuracy of these 
panel method code solutions (see App. 4 and App. 5). 

The calculation of wall interference in an open jet using panel method code ANTARES 
is similar to calculations described in the previous chapter. Coefficients describing the wall 
boundary condition on the four test section walls are assigned such that the open jet 
boundary condition is imposed. The following table lists values of these coefficients for the 
open jet boundary condition : 


Table 3 : Wall Boundary Condition Coefficients (Open Jet) 


Wall No. 

Panel Index Range 

Cl(») 

C2 (*) 

c 3 (i) 

c A (i) 

1 

1 < i < 1200 

0 

1 

0 

0 

2 

1201 < t < 2400 

0 

1 

0 

0 

3 

2401 < i < 3600 

0 

1 



4 

3601 < i < 4800 

0 

1 
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An open jet with square cross-section is selected for the study of the accuracy of panel 
method code ANTARES. An open jet width and height of 10 [ft] is selected. The open 
jet boundaries are represented by using 4800 panels (60 panels in the streamwise direction 
and 80 panels in the circumferential direction, see Fig. 4). The paneled test section length 
is equal to eight tunnel widths. Compressibility effects are investigated by computing wall 
interference for two different Mach numbers (0.00 and 0.90). 

A single point doublet or line doublet is placed inside of the tunnel to represent 
blockage or lifting effects of the wind tunnel model (see Fig. 3). Panel method code 
ANTARES and the Method of Images are used to compute the wall interference flow field 
on the tunnel centerline (see Line 1 in Fig. 3). 

Figures 11a, lib, and 11c depict the result of the wall interference calculation for 
the selected unit point doublet representing test article blockage in an open jet. The 
dimensionless wall interference perturbation velocity components are depicted. The overall 
agreement between panel method code solution (abbreviated P.M.) and Method of Images 
solution (abbreviated M.I.) for both Mach numbers is excellent. Differences between panel 
method code solution and Method of Images solution of the axial interference velocity 
component are smaller than for a corresponding calculation using the closed wall boundary 
condition (compare Fig. 11a with Fig. 5a) 

Figures 12a, 12b, and 12c depict the result of the wall interference calculation for a 
unit line doublet representing test article lifting effects in an open jet. The line doublet 
orientation angle is selected to be 0.0° . The dimensionless wall interference perturbation 
velocity components are depicted. Again, the agreement between panel method code 
solution (abbreviated P.M.) and Method of Images solution (abbreviated M.I.) for both 
Mach numbers is excellent. 

2.6 Perforated Wall Boundary Condition 

Classical solutions of blockage and lift interference effects in a two-dimensional perfo- 
rated wall wind tunnel may be used to check the application of perforated wall boundary 
conditions in panel method code ANTARES. The classical solution of the perturbation 
potential and the blockage factor caused by a two-dimensional point doublet on the cen- 
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terline of a two-dimensional perforated wall wind tunnel is given by Baldwin et al. [6]. 
Similarly, the classical solution of the lift interference caused by a two-dimensional vortex 
on the centerline of a two-dimensional perforated wall wind tunnel is given by Pindzola 
and Lo [7]. 

In the first part of this chapter, classical solutions are described in greater detail. 
Afterwards classical and panel method code solutions for blockage and lift interference are 
compared and discussed. 

Classical Blockage Factor : Baldwin et al. [6] use the Fourier transform technique 
to compute the perturbation potential <p w due to tunnel walls (<p* in their nomenclature). 
The blockage factor on the tunnel centerline may easily be computed by differentiating 
Eq. (19) of Ref. [6] with respect to the streamwise coordinate. It is assumed that the walls 
of the two-dimensional tunnel are “ideal porous” (K = 0) and that wall interference is 
computed on the tunnel centerline ( y = 0). Then. Eq. (19) of Ref. [6] becomes : 


<Pw(x) 


— m e 
2 7T 0 h 



h{q) dq + 



h{q) d q 


(24 a) 


where 

j = _L i 

R ( cosh(q) ] 2 + [ (P/R) ■ sinh(q) ] 2 

T _l_ . [ l - (P/Rf ] + [ l + tf/R) 2 ] - e~ 2 « 

2 [ cosh(q) ] 2 + [ (P/R) • sinh(q) ] 2 \Ph ) \ ' 

Parameter m e in Eq. (24a) may be expressed as a function of the point doublet 
strength a by using Eq. (2.43) from Ref. [2] or Eq. (11) from Ref. [6]. Then, we get : 



m e 


' ^oo 

~ — 



(25) 


The perturbation potential due to tunnel walls (Eq. (24a)) may be differentiated with 
respect to the streamwise coordinate x. Combining the corresponding result with Eq. (25) 
and with the definition of the blockage factor we get : 
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where 

d _ h = JZA. . — r • sin ( (266) 

d x Rh [ cosh(q) f + [ (0/R) • sinh(q) ] \ 0 h / 

dh _ _j_ . [ i - W-R) 2 ] + 1 1 + W*> 2 1 ■ e ~ 2 ’ . cos ( -UL) (26c) 

d x 2 0 h [ cosh(q) f + [ (0/R) • sinh(q) ] 2 \ 0 h / 

Baldwin et al. [6] give a solution for the blockage factor at the location of the two- 
dimensional point doublet (see Eq. (22) in Ref. [6]). This solution may be used to check 
the above result. Using the above equations, we get for blockage factor e at the location 
of the point doublet (x = 0) : 


e (0) 


1 

^oo 
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2 7T 3 2 h 
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2 7T 0 h 
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where 
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1 [ 1 - (0/R) 2 ] + [ 1 + (0/R) 2 ] ■ e" 2? 

2 2 


0 h [ cosh(q) f + [ (0/R) ■ sinh(q) ] 

We also know that : 
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2 2 

= (cosh(q) - (0 / R) 2 ■ sinh(q) ) ■ e~ q ■ q 

Combining Eqs. (27a), (27b), and (28), we get the blockage factor at the location of the 
point doublet : 


e(0) = 


— m. 


Uoo 2 n 0 2 h 2 


j; 


[ cosh(q ) - ( 0/Rf ~ sinh(q) ] • e q _ • g ^ ^ 


[ cos%) ] 2 + [ (/?/#) • «»»/»(«) ] 


Equation (29) is equal to Eq. (22) of Ref. [6]. Thus, Eqs. (26a), (26b), (26c) reduce to 
Eq. (22) of Ref. [6] if the blockage factor is computed at the location of the two-dimensional 

point doublet. 

Classical Lift Interference : Pindzola and Lo [7] also use the Fourier transform 
technique to compute the lift interference caused by a two-dimensional vortex in a two- 
dimensional perforated wall wind tunnel. They express lift interference in terms of upwash 
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velocity w w (w in their nomenclature). Using Eq. (3.65) of Ref. [7], the dimensionless 
upwash perturbation velocity component caused by a two-dimensional vortex on the cen- 
terline (z = 0) of a two-dimensional perforated wall wind tunnel may be written as : 
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W 


Uoo 


i - r 

Uoo 2 7T h 



r°° 

I I 4 d q + 

/ hdq 
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where 


q x 

Th 


I [ sinh(q) — ( 0 / R ) 2 ■ cosh(q) ■ e~ q 

[ sinh(q) ] 2 + [ ( 0 j R) • cosh(q) ] 2 
Comparison of Wall Interference : Boundary conditions in panel method code 
ANTARES may be modified so that the flow field of a two-dimensional tunnel is computed. 
It is only required to impose the ” closed wall” boundary condition on the tunnel sidewalls 
(Wall 1 and Wall 2 in Fig. 1) and the ’’perforated wall” boundary condition on the tunnel 
floor and ceiling (Wall 3 and Wall 4 in Fig. 1) of a tunnel with square cross-section. Table 4 
lists wall boundary condition coefficients selected for the simulation of a two-dimensional 
perforated wall wind tunnel (see Eq. (1)). 


q x 
~0h 
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Table 4 : Wall Boundary Condition Coefficients (2D Perforated Wall) 


Wall No. 

Panel Index Range 

Cl (0 

C2 (»') 

c 3 (0 

C 4(0 

1 

1 < i < 1200 

0 

0 

1 

0 

2 

1201 < i < 2400 

0 

0 

1 

0 

3 

2401 < i < 3600 

0 

1 

1 /R 

0 

4 



1 




The potential and derivatives of the three-dimensional point doublet in panel code 
ANTARES have to be replaced by corresponding equations of a two-dimensional point 
doublet (see Ref. [8], p.67/68). It is also required to replace the potential and derivatives 
of the three-dimensional line doublet by corresponding equations of a two-dimensional 
vortex (see Ref. [8], p.70/71). 

A square tunnel (10 [ft] x 10 [/<], see Fig. 3 and Fig. 4) is selected for the comparison 
between classical solutions (Eqs. (26a), (26b), (26c) for blockage and Eqs. (30a), (30b), 
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(30c) for lift interference) and corresponding panel method code solutions computed by 
ANTARES. The restriction parameter R of the perforated floor and ceiling of the square 
tunnel is set to 1.14 (restriction parameter of NASA Ames lift TWT, see Ref. [2], p.3-25). 
Two different Mach numbers (0.00 and 0.90) are selected to investigate compressibility 

effects on the wall interference calculation. 

At first, blockage corrections are investigated. A two-dimensional unit point doublet 
pointing upstream is placed on the tunnel centerline at the coordinate system origin (<r 
l x _ y _ z _ o). Results of the panel method code calculation are compared in Fig. 13 
with the corresponding exact solution (Eqs. (26a), (26b), (26c)). The agreement between 
panel method code solution and exact solution is excellent for both Mach numbers. 

In the next step, upwash corrections are investigated. A two-dimensional unit vortex 
is placed on the tunnel centerline at the coordinate system origin (T = 1, x = y = 2 = 0). 
Results of the panel method code calculation are compared in Fig. 14 with the correspond- 
ing exact solution (Eqs. (30a), (30b), (30c)). Again, the agreement between panel method 
solution and exact solution is excellent for both Mach numbers. 
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CHAPTER 3 


CONCLUSION AND REMARKS 

Panel method code ANTARES was developed to determine wall interference effects 
of point doublets and line doublets in a wind tunnel with rectangular cross-section. The 
code is based on Keller's algorithm [3] that allows the user to describe up to six different 
types of wall boundary conditions. 

At the present time ANTARES can assign the closed wall, open jet, or perforated 
wall boundary condition to the centroid of each wall panel. Blockage effects of the test 
article are represented by using point doublets. Test article lifting effects are represented 
by using line doublets. Compressibility effects are taken into account by applying the 
Prandtl-Glauert transformation to the test section geometry, the point doublet or line 
doublet coordinates, and to the singularity strength. 

The accuracy of the panel method code calculation for a given panel representation 
of the test section geometry was investigated by using available classical solutions. Panel 
method code solutions of the closed wall and open jet boundary condition were successfully 
compared with corresponding Method of Images solutions. Panel method code solutions 
of the perforated wall boundary condition were successfully compared with corresponding 
Fourier transform solutions. 

Studies have shown that the accuracy of a wall interference calculation depends on the 
number of wall panels used to represent the test section walls. The accuracy also depends 
on the free— stream Mach number that is used to transform the test section geometry. 
For a fullspan model test section configuration, a total number of 4800 panels (60 in the 
streamwise direction and 80 in the circumferential direction) seems to provide a good 
representation of the test section geometry. For a semispan model test, a total number 
of 7200 panels (60 in the streamwise direction and 120 in the circumferential direction) 
should be selected. In both cases it is recommended to use cosine panel spacing in the 
streamwise direction and equal panel spacing in the circumferential direction. Accuray 
improvements were achieved by increasing the number of panels in the circumferential 
direction. The length of the paneled test section should be about eight tunnel widths. 
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ANTARES computes sufficiently accurate corrections up to a Mach number of 0.80 if the 
closed wall boundary condition is applied to the selected panel representation of the wind 
tunnel walls. The Mach number limit seems to be close to 0.90 if the open jet or perforated 
wall boundary condition is applied. 

Panel method code ANTARES does not model the effect of shocks that may occur 
during the test of a wind tunnel model at high subsonic Mach numbers. Mach number 
limits for the closed wall, open jet, or perforated wall boundary condition are therefore 
only valid as long as no shocks touch the test section walls during the wind tunnel test. 

Derivatives required for the application of the ideal slotted wall boundary condition 
and of the slotted wall boundary condition that includes viscous effects in slots have al- 
ready been included in ANTARES (see Table 1 on p.3 and Appendices 1 and 2). An 
implementation of these boundary conditions in ANTARES will be considered in a future 
release of the software. 
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APPENDIX 1 


POINT DOUBLET POTENTIAL AND DERIVATIVES 


Panel method code ANTARES uses point doublets pointing upstream to describe 
blockage effects of a test article in a rectangular wind tunnel. ANTARES describes wind 
tunnel wall boundary conditions by applying Keller’s algorithm (see Ref. [3]). Therefore, 
it is necessary to specify the point doublet potential and its derivatives in the tunnel 
coordinate system x,y,z . 

A point doublet is located at point x 0 ,y 0 ,z 0 in the tunnel coordinate system. The 
point doublet is pointing upstream, i.e. in the negative x-direction (see Fig. 1). Then, the 
potential of the point doublet may be written as [8]: 


<f>T 


X — Xc 


4 7 r 


(x - x 0 ) 2 + (y - y 0 ) 2 + (z - Zo) 2 


3/2 


( 1 . 1 ) 


First and second order derivatives may easily be obtained by differentiating Eq. 
For the first order derivatives we get : 


( 1 . 1 ). 


d qi> m 

d x 


a 

4 7 r 


( y-y° ) 2 + (z - z 0 ) 2 - 2 • (x — x 0 ) 2 


(x -x, ) 2 + (y - y e ) 2 + (z - z o y 


I 5/2 


(1.2a) 


d 4> m 

(-3)-a 

(x - Xo) • (y - y c ) 


d y 

4 * 7T 

(x - Xo) 2 + (y - yo ) 2 + (z - z 0 ) 2 

5/2 

d 4> m 

(-3) -a 

(x - x 0 ) • (z - Zo) 


d z 

4 • 7 r 

(x - Xo) 2 + (y - y Q ) 2 + (z - Zo) 2 

5/2 


Similarly, for the second order derivatives we get : 


(1.26) 


(1.2c) 


d 2 

d x d y 


d 2 <ft m 

d x d z 


(_ 3 ) . a (y-y°)- (y~y °) 2 + (---°) 2 - 4-(x-x 0 ) 2 


4 * 7T 


(x - Xo) 2 + (y - yo) 2 + {z - z o y 


U/2 


(1.3a) 


(_ 3 ) . a (z-zo)- (y~y o) 2 + (z-z 0 ) 2 - 4 • (x — x 0 ) 2 


4 * 7T 


(x - Xo) 2 + (y - y Q ) 2 + (z - z 0 ) 2 


1 7/2 


(1.36) 
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Derivatives normal to the walls of the rectangular wind tunnel are required as panel 
method code ANTARES uses Keller’s algorithm [3] to model different types of wall bound- 
ary conditions. In general, the normal derivative operator may be written as : 


d n 


= n o V = 




The scalar product defined in Eq. (1.4) may also be written as : 

d d __ d 


d n 


= N t 


d x 


+ N v 


d y 


+ N z 


d z 


(1.4) 


(1.5) 


Now it is possible to determine the normal derivatives if the outward normal vector 
n on each wind tunnel wall is identified. Using Fig. 1 and Eq. (1.5), we get : 


Table 5 : Normal Derivatives on the Wind Tunnel Wall (Point Doublet) 


Wall No. 

n 

d 

dn 

d 4> m 

dn 

d 2 4>m 

d x d n 

1 


W 


d 

d y 

d <j>m 
d y 

Eq.(1.2b) 

d 2 <j> m 

d x d y 

Eq.(1.3a) 

2 

t 

0 N 
-1 
0 , 


d 

d y 

d 4> m 
d y 

Eq.(1.2b) 

d 2 4> m 

d x d y 

Eq.(1.3a) 

3 


(!) 


d 

d z 

d <f>m 

d z 

Eq.(1.2c) 

d 2 <t> m 

d x d z 

Eq.(1.3b) 

4 

( 

( 0 ' 

0 

) 

d 

d z 

d <f>m 
d z 

Eq.(1.2c) 

d 2 4>m 
d x d z 

Eq.(1.3b) 
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APPENDIX 2 


LINE DOUBLET POTENTIAL AND DERIVATIVES 

Panel method code ANTARES uses line doublets to model lifting effects of a test 
article in a rectangular wind tunnel. Wind tunnel wall boundary conditions axe described 
by applying Keller’s algorithm [3]. Therefore, it is necessary to specify the line doublet 
potential and its derivatives in the tunnel coordinate system x,y,z . 

Figure 15a shows the location of a line doublet in a rectangular wind tunnel. The line 
doublet is oriented such that the corresponding lift force points in the £ direction of the 
line doublet fixed coordinate system £,?7,£ . Line doublet orientation angle 0 is used to 
describe the line doublet orientation in the tunnel coordinate system x,y,x . 

The calculation of the line doublet potential and its derivatives in the tunnel coordi- 
nate system x,y,z may be done in three steps: 

(1) A coordinate system transformation from the tunnel coordinate system x,y,x to 
the line doublet fixed coordinate system £,77,£ is introduced which will greatly simplify the 
calculation of the line doublet potential and its derivatives. 

(2) The line doublet potential and its derivatives are determined in the line doublet 
fixed coordinate system £,t/,£ . 

(3) The line doublet potential and its derivatives are transformed from the line doublet 
fixed coordinate system £,77,£ to the tunnel coordinate system x,y,z . 

Coordinate System Transformation 

In general, it is easier to compute the line doublet potential and its derivatives in 
the line doublet fixed coordinate system £,77, £ . Therefore, the following transformation 
may be used to transform coordinates from the tunnel coordinate system x,y,z to the line 
doublet fixed coordinate system £,?7,£ . Thus, we get for a flow field point x,y,z in the 
tunnel coordinate system : 


£(*) 

= X 



(2.1a) 

7(y>2,0) 

= y ■ 

cos © — z • 

■ sin © 

(2.16) 

C(y,x,0) 

= y ■ 

sm Q + z 

• cos 0 

(2.1c) 


Similary, starting point coordinates x 0 ,y 0 ,2 0 of the line doublet have to be transformed to 
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the line doublet fixed coordinate system (,vX . We get : 


(2.2a) 

(2.26) 

(2.2c) 


Co — Xo 

Tj 0 = y 0 • cos 0 — z 0 ■ sin 0 

£ 0 = y 0 ■ sin 0 + z 0 ■ cos 0 

Potential and Derivatives in Line Doublet Fixed Coordinate System 

The line doublet potential in the line doublet fixed coordinate system £,t?,C may be 
written as [1]: 


0m 


a 

4 7T 


C-Co 

(V - ^o ) 2 + (C - Co ) 2 


• A 


(2.3 a) 


where 


A = 1 + 


C ~ Co 

(C - Co ) 2 + {v - Vo) 2 + (C - Co ) 2 


I 1/2 


(2.36) 


Derivatives of the line doublet potential in the line doublet fixed coordinate system 
may now be obtained by differentiating Eq. (2.3a). For the first order derivatives we 

get : 


d 4> m 

di 


c r 

4 7 r 


C ~~ Cq 

(C - Co ) 2 + (v~ Vo) 2 + (C - Co ) 2 



(2.4) 


d <t>rn 

d r] 


where 


- <y_ _ iv - Vo) • (C ~ Cq) 

4 7 r (r? - Vo) 2 + (C - Co ) 2 


2_B 

(v - Vo) 2 + (C - Co ) 2 


+ c 


(2.5 a) 


B = 1 + 


£ ~ 

(C ~ Co ) 2 + (v ~ Vo) 2 + (C — Co ) 2 


1 1/2 


(2.56) 
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Potential and Derivatives in Tunnel Coordinate System 

The line doublet potential and its derivatives in the tunnel coordinate system may eas- 
ily be obtained by reversing the coordinate system transformation defined by Eqs. (2.1a) to 
(2.1c). The line doublet potential is a scalar function and depends only on the coordinates 
of the selected flow field point. Thus, combining Eqs. (2.1a) to (2.1c) with Eqs. (2.3a), 
(2.3b), we get for the potential : 

<t> m (x,y,z) = <t> m (Z(x),r}(y,z,Q), ((y,z,Q)) (2.9) 
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Derivatives in the tunnel coordinate system may be obtained by reversing the coordi- 
nate system rotation (see Fig. 15b). Then, we get for the first order derivatives of the line 
doublet potential : 


d <p 


m 


d <f>rn 
d x 

d 4>m 


d y 


d rj 


d <t>m 

d£ 

cos 0 + 


d <$>, 


0<f> r 


0 4> r 


( —sin 0 ) + 


d C 

6 <t>; 


• sin 0 


cos 0 


(2.10a) 

(2.106) 

(2.10c) 


dz drj v ; dC 

The higher order derivatives may be obtained in a similar manner. We know from 
Eq. (2.1a) : 




Thus, we get 


d 2 4>r 


d x d y 


d 2 4>, 


= X =* 
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d£ 

-<2 

III 

(2.11a) 

d 2 <j>m 

d 

d 4>m 

(2.116) 

d x d y 

d£ 

d y 

d 2 <Pm _ 

d 

d 4>m 

(2.11c) 

d x d z 

i (2.11b), we get 

d£ 

d z 

d 2 4>m 

d £ d t) 

i (2.11c), we get 

^ , d 2 <j>m ^ 

■ cos 0 + . _ . • Sin 0 

d£d C 

(2.12a) 

d 2 4>m f r\\ > 

di a n ' ( sm0) + 

d(d c c “® 

(2.126) 


d x d z 

Derivatives normal to the walls of the rectangular wind tunnel are required as panel 
method code ANTARES uses Keller’s algorithm [3] to model different types of wall bound- 
ary conditions. From App. 1 we know that the normal derivative operator is equal to : 


d n 


= N x 


d 


d x 


+ N v 


d y 


+ N z 


d 


d z 


(2.13) 


Now it is possible to determine the normal derivatives if the outward normal vector 
n on each wind tunnel wall is identified. Using Fig. 1 and Eq. (2.13), we get : 
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Table 6 : Normal Derivatives on the Wind Tunnel Wall (Line Doublet) 


Wall No. 

n 

d 

d <Pm 

d 2 (p m 

d n 

d n 

d x d n 
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(!) 

d 

d y 

d 4>m 

d y 

Eq.(2.10b) 

d 2 <p m 

d x d y 

Eq.(2.12a) 
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■ 

d 

d y 

9 0 m 
d y 

Eq.(2.10b) 

d 2 <f> m 

d x d y 

Eq.(2.12a) 

3 


(!) 

d 

d z 

d <pm 

d z 

Eq.(2.10c) 

d 2 <t>m 

d x d z 

Eq.(2.12b) 

4 


1 

d 

d z 

d 4> m 

dz 

Eq.(2.10c) 

d 2 4>m 

d x d z 

Eq.(2.12b) 


The calculation of the line doublet potential and its derivatives in the tu nn el coordi- 
nate system is complex as the line doublet orientation angle has to be taken into account. 
A coordinate system transformation from the tunnel coordinate system to a line doublet 
fixed coordinate system is required. Then, the potential and its derivatives are computed in 
the line doublet fixed coordinate system. Finally, the potential and its derivatives are com- 
puted in the tunnel coordinate system by reversing the coordinate system transformation. 
The following flow chart summarizes different solution steps : 


41 























Selection of flow field point (x,y,z) in tunnel coordinate system. 
Selection of line doublet (x o ,j/o,2 o ,0) in tunnel coordinate system. 

•n 

Transformation from tunnel coordinate system to line 
doublet fixed coordinate system (Eqs. (2.1a) to (2.2c)) 

Calculation of potential and its derivatives in line 
doublet fixed coordinate system (Eqs. (2.3a) to (2.8)) 

Calculation of potential and its derivatives in 
tunnel coordinate system (Eqs. (2.9) to (2.12b)) 
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APPENDIX 3 


SINGULARITY STRENGTH TRANSFORMATION 


In general, the Prandtl-Glauert transformation has to be applied if the subsonic wind 
tunnel and wall interference flow field of a singularity located in a wind tunnel is calculated. 
This transformation requires a stretching of the wind tunnel coordinates from an original 
coordinate system ( x,y,z ) to a corresponding transformed coordinate system (x' ,y',z'). 
The transformation of the coordinates, of the free— stream velocity, and of the perturbation 
potential is given by the following equations (see Ref. [4], p.217-223 and p.348-351) : 


X 

= X 

(3.1a) 

y' 

= y-P 

(3.16) 

t 

Z 

= z-p 

(3.1c) 

<4 

— W CO 

(3.1d) 

4>' 

= 4>-P 2 

(3-le) 

P 

= V i - ML 

(3-1/) 


where Ucc is the free-stream velocity, <f> is the perturbation velocity potential, and is 
the free-stream Mach number of the subsonic flow field. After the application of the trans- 
formation, the flow field solution is found in the transformed coordinate system by solving 
the Laplace equation. Finally, a back-transformation has to be applied to the solution of 
the Laplace equation in the transformed coordinate system to obtain the solution of the 
subsonic potential equation in the original coordinate system. 

Singularities are used to represent a test article during a wind tunnel wall interference 
calculation. Point sources or point doublets represent blockage effects of the test article. 
Vortices or line doublets represent lifting effects of the test article. The geometry of the test 
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article is related to the strength of these singularities. The Prandtl-Glauert transformation 
changes the test article geometry. Thus, the Prandtl-Glauert transformation also has to 
change the strength of singularities representing the test article. The transformation of 
the singularity strength has to be a function of scale factor /?, i.e. 


a' = <r-F(p) (3.1 g) 


The particular form of the transformation function F(/3 ) for each type of singularity 
is rigorously derived in the following parts of this section. 

Two-Dimensional Point Source 

A two-dimensional point source may be used to describe the blockage effect of a 
halfbody in a two-dimensional flow field. The strength cr of this point source in the 
original coordinate system ( x,z ) may be written as (see Ref. [9], p.62) : 


cr — Uqo ■ D (3-2) 


where is the free-stream velocity and D is the width of the two-dimensional halfbody. 
The singularity strength in the transformed coordinate system (x' ,z') may be written as : 

cr' = • £>' (3.3) 

Applying the Prandtl-Glauert transformation, i.e. Eq. (3.1c) and (3.1d), to the width 
of the two-dimensional halfbody and to the free-stream velocity, we get : 

u'oo ’D' = Uoo •£>•/? (3.4) 

Finally, combining Eqs. (3.2), (3.3), and (3.4), we get : 

cr' = a ■ (3 (3.5a) 

where a ' is the transformed strength and cr is the original strength of a two-dimensional 
point source. Thus, for a two-dimensional point source, the transformation function F(J3) 
as defined in Eq. (3.1g) is given as : 
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F(0) = P (3.56) 


It is possible to verify the singularity strength transformation given in Eq. (3.5b). 
Therefore, a Method of Images solution for the blockage factor caused by a 2D point 
source at a high subsonic Mach number was developed using derivatives of the singularity 
potential in combination with Eqs. (3.5a), (12a), (12c), and (23a). This solution for the 
blockage factor showed good agreement with a classical solution given by Eq. (2.56) in 
Ref. [2], 

Two-Dimensional Point Doublet 

A two-dimensional point doublet may be used to represent blockage effects of a cylin- 
der in a two-dimensional flow field. The strength a of this point doublet may be written 
as (see Ref. [2], Eq. (2.43) or Ref. [8], Eq. (3.96)) : 


<7 = Moo • 2 • 7T • r 2 c (3.6) 


where is the free-stream velocity and r c is the radius of the cylinder. The product 
x ' r c i n Eq. (3.6) may be interpreted as the cross-sectional area of the cylinder. This area 
can also be written in integral form as : 


7 r ■ r! 


r+ r c 

= / z(x) d x (3.7) 

^ ~ T C 

where z(x ) is the width of the cylinder for a given streamwise location x . Combining Eqs. 
(3.6) and (3.7), the singularity strength in the original coordinate system ( x,z ) may be 
written as : 


— * 2 


/■+ r c 

J—r r 


z(x) d x 


(3.8) 


Thus, the singularity strength in the transformed coordinate system (x 1 ,z') may be 
written as : 

r+ r c 

= u 'oo- 2 J_, *V) d *' (3.9) 
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Applying the Prandtl-Glauert transformation, i.e. Eqs. (3.1a) and (3.1c), to the 
streamwise coordinate x and the width z(x) of the cylinder we get : 


a = 


d x =>■ r' c = 

r c 

(3.10a) 

'(*') = z{x)-0 


(3.106) 

si), and (3.10b), we get : 



r+r c 



Uoc ■ 2 • / z{x ) d x 

J-r, 

■3 

(3.11) 


Combining Eqs. (3.8) and (3.11), we get : 

a' = <7 • 0 


(3.12a) 


where a f is the transformed strength and a is the original strength of a two-dimensional 
point doublet. Thus, for a two-dimensional point doublet, the transformation function 
F(I3) as defined in Eq. (3.1g) is given as : 



It is possible to verify the singularity strength transformation given in Eq. (3.12b). 
Therefore, a Method of Images solution for the blockage factor caused by a 2D point dou- 
blet at a high subsonic Mach number was constructed using derivatives of the singularity 
potential in combination with Eqs. (3.12a), (12a), (12c), and (23a). This solution for the 
blockage factor showed good agreement with a classical solution given by Eq. (2.45) in 
Ref. [2]. 

Three-Dimensional Point Source 

A three-dimensional point source and point sink may be used to represent solid volume 
blockage effects of a fuselage in a three-dimensional flow field. The strength a of a three- 
dimensional point source may be written as (see Ref. [9], p.65) : 
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a 


(3.13) 


= Uoo • 7T • 

where is the downstream radius of the three-dimensional body of revolution that 
corresponds to a single source in a flow field and is the free-stream velocity. The 
downstream radius of the three-dimensional body of revolution may be written as : 

rlo = Vi*? + z{*) 2 (3.14) 

where y(x) and z(x) describe the cross-sectional shape of the body of revolution far down- 
stream of the source. Combining Eqs. (3.13), (3.14) we get the strength a of a three- 
dimensional point source in the original coordinate system ( x,y,z ) : 

a = Uoo • 7r • y{x) 2 + z(x) 2 J (3.15) 

The singularity strength in the transformed coordinate system (x',y',z') may be writ- 
ten as : 

a ' = u'oo • * • y \ x ') 2 + z '{ x ') 2 (3.16) 

Applying the Prandtl-Glauert transformation, i.e. Eqs. (3.1a), (3.1b), (3.1c), we get : 


II 

(3.17a) 

z'(x') = z(x) ■ (3 

Combining Eqs. (3. Id), (3.16), (3.17a), (3.17b), we get : 

(3.175) 

a' = Uoo-tt- y{x) 2 + z(x) 2 ■ (3 2 

(3.18) 

Finally, combining Eqs. (3.15) and (3.18), we get : 


a' = a ■ (3 2 

(3.19a) 

where o' is the transformed strength and a is the original strength of a three-dimensional 
point source. Thus, for a three-dimensional point source, the transformation function 
F(/9) as defined in Eq. (3.1g) is given as : 
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It is possible to verify the singularity strength transformation given in Eq. (3.19b). 
Therefore, a Method of Images solution for the blockage factor caused by a 3D point 
source at a high subsonic Mach number was constructed using derivatives of the singularity 
potential in combination with Eqs. (3.19a), (12a), (12b), (12c), and (23a). This solution 
for the blockage factor showed good agreement with a classical solution given by Eq. (2.72) 
in Ref. [2]. 


Three-Dimensional Point Doublet 

A three-dimensional point doublet may be used to represent solid volume blockage 
effects of a small fuselage in a three-dimensional flow field. The strength cr of a three- 
dimensional point doublet may be written as (see Ref. [8], p.80) : 



where r 0 is the radius of the sphere and u 0 0 is the free-stream velocity. The product 2-n -r\ 
in Eq. (3.20) may be interpreted as a multiple of the volume of a sphere. Therefore, we 
can write : 



where y{x),z{x) describes the cross-sectional shape of the sphere as a function of the 
streamwise coordinate. Combining Eqs. (3.20) and (3.21), we get for the singularity 
strength in the original coordinate system (x,y,z) : 



(3.22) 


The singularity strength in the transformed coordinate system (x' ,y' ,z') may be writ- 
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ten as : 


, , 3 f+< 



* = ■ T ■ y_ r ; - 

y\x) 2 + z'(x') 2 j d x 

(3.23) 

Applying the Prandtl-Glauert transformation, 

i.e. Eqs. (3.1a), (3.1b), (3.1c), we get : 

d x’ = d x =* 

t 

> T = TV 

' O '0 

(3.24a) 

yV) = 

y(x) ■ (3 

(3.246) 

*V) = 

z{x) ■ 0 

(3.24c) 


Combining Eqs. (3.1d), (3.23), (3.24a), (3.24b), (3.24c), we get : 


II 

£ 

8 

o| CO 

r+r 0 
/ 7r ‘ 

y(x ) 2 + z (x ) 2 

1 

H 

^3 

1 

L 

0 

1 

L -* 

1 J 


Finally, combining Eqs. (3.22), (3.25), we get : 

= CT ‘/? 2 (3.26a) 

where a is the transformed strength and a is the original strength of a three-dimensional 
point doublet. Thus, for a three-dimensional point doublet, the transformation function 
F((3) as defined in Eq. (3.1g) is given as : 


F(0) = (3 2 (3.266) 


It is possible to verify the singularity strength transformation given in Eq. (3.26b). 
Therefore, a Method of Images solution for the blockage factor caused by a 3D point dou- 
blet at a high subsonic Mach number was constructed using derivatives of the singularity 
potential in combination with Eqs. (3.26a), (12a), (12b), (12c), and (23a). This solution 
for the blockage factor showed good agreement with a classical solution given by Eq. (2.64) 
in Ref. [2]. 
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Two-Dimensional Vortex 

A two-dimensional vortex may be used to represent lifting effects of an airfoil in a 
two-dimensional flow field. The strength a of a two-dimensional vortex may be derived by 
using the Kutta-Joukowski formula. The Kutta-Joukowski formula may be written as : 

L /As = PooUoo-T (3.27) 


where L/As is the lift force per unit wing span, poo is the free-stream density, Uoo is the 
free — stream velocity, and T is the circulation. Thus, in case of a two-dimensional vortex, 
we may define the strength as : 


= T (3.28) 

Poo ■ Uqo ■ As 


The lift force per unit wing span, i.e. L/As, may also be written by using the lift 
coefficient definition. Then, we get : 

L/As = -£=- ■ • cl • e (3.29) 

where c is the chord of the airfoil. Combining Eqs. (3.28), (3.29) we get : 


a 


tioo 

T~ 


CL - C 


(3.30) 


The lift coefficient c L may be rewritten by considering the pressure coefficient differ- 
ence Ac p between upper and lower surface of an airfoil as a function of the streamwise 
coordinate x. Using Eq. (2-10) from Ref. [4], we get : 


cl = 



(3.31) 


Combining Eq. (3.30) and (3.31), we get for the vortex strength in the original coordinate 
system (x,y,z) : 

a = • f° Acp(x) d x (3.32) 

^ J 0 


50 



The singularity strength in the transformed coordinate system (x',y r ,z') may be writ- 
ten as : 



The stream wise coordinate x and chord c are transformed by using Eq. (3.1a) : 


dx’ = dx => c' = c (3.34a) 

The transformed pressure coefficient difference may be written as a function of the 
original pressure coefficient difference. Assuming small perturbations, we know for the 
pressure coefficient in the original coordinate system : 


c P (x ) = 


-2 d<f> 


U r 


•-oo ^ X 

Similarly, we know for the pressure coefficient in the transformed coordinate system : 

-2 dd>' 


(3.346) 


c'(x') = 


u 'ao 9 X' 


(3.34c) 


Using Eq. (3.1a), (3. Id), (3.1e) and applying the chain rule to Eq. (3.34c), we get : 


c' P (x') 


-2 d x d ( d> • (3 2 ) 


U oo d x' d x 


-2 d 4> 


Ur 


d x 


■P 2 


(3.34d) 


Combining Eqs. (3.34b) and (3.34d), we get : 


C ' P ( X ') = c p (x)-0 2 (3.34e) 

Equation (3.34e) agrees with Eq. (4-22) from Ref. [4]. Thus, we may write the pressure 
coefficient difference in the transformed coordinate system as r 


A C ;(x') = A c p (x)-0 2 


(3.34/) 


Combining Eqs. (3.1d), (3.33), (3.34a), and (3.34f), we get : 



(3.35) 
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Finally, combining Eqs. (3.32), (3.35), we get : 

a' — cr ■ (3 2 ( 3.36a) 

where cr' is the transformed strength and cr is the original strength of a two-dimensional 
vortex. Thus, for a two-dimensional vortex, the transformation function F(f3) as defined 
in Eq. (3.1g) is given as : 


F{0) = /? 2 (3.366) 


It is possible to verify the singularity strength transformation given in Eq. (3.36b). 
Therefore, a Method of Images solution for the upwash velocity caused by a 2D vortex at a 
high subsonic Mach number was constructed using derivatives of the singularity potential 
in combination with Eqs. (3.36a), (12a), (12c), and (23c). This solution for the upwash 
velocity showed good agreement with a classical solution given by Eq. (2.13) in Ref. (2). 

Three-Dimensional Line Doublet 

A three-dimensional line doublet may be used to represent lifting effects of a wing 
in a three-dimensional flow field. The strength cr of a three-dimensional line doublet is 
defined as lift force divided by density and free-stream velocity : 



L 

<0 

<1 

u 

II 

(3.37) 

<J — 

Poo ’ ^OO 




where is the free-stream velocity, poo is the density, T is the circulation, and As is 
the span of a wing span element. The lift force L, may also be written by using the lift 
coefficient definition. Then, we get : 

L = (3.38) 

where c is the chord of the airfoil. Combining Eqs. (3.37), (3.38) we get . 

a — — * ci * c * A s (3.39a) 
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The lift coefficient may be expressed as a function of the pressure coefficient difference 
on the surface of the wing span element. Using Eq. (3-54a) from Ref. [4], we get : 

Y /*c rAs/2 

CL = c-A S Jo J-As/2 ACp ^ ,y ) dx dy (3.396) 

Combining Eq. (3.39a) and (3.39b), we get for the line doublet strength in the original 
coordinate system (x,y,z) : 


■ „ . 2 =. . r f As/ 2 

2 JO J-As/2 

A c p (x,y) d x d y 

(3.40) 

The singularity strength in the transformed coordinate system (x',y' r 

z') may be writ- 

ten as : 



✓ . . r' f As ' r - 

2 JO J-As'/2 

^ c p( x ^y , ) d x' d y' 

(3.41) 

The y coordinate and wing span As are transformed by using Eq. (3.1b). 

Then, we get : 

d y' = d y-p =» 

A s' = As • P 

(3.42) 


Combining Eqs. (3.1d), (3.34a), (3.34f), (3.41), and (3.42), we get : 

, Uoo f C f As / 2 1 

a “ — 7 0 j_ As /2 &c p (x,y) d x dy ■ p 3 (3.43) 

Finally, combining Eqs. (3.40), (3.43), we get : 

= cr ■ (3 3 (3.44a) 

where a' is the transformed strength and a is the original strength of a three-dimensional 
line doublet as defined in Eq. (3.3 i ). Thus, for a three-dimensional line doublet strength 
as defined in Eq. (3.37), the transformation function F(j3) as defined in Eq. (3.1g) is given 
as : 
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It is possible to verify the singularity strength transformation given in Eq. (3.44b). 
Therefore, a Method of Images solution for the upwash velocity caused by a 3D line dou- 
blet at a high subsonic Mach number was constructed using derivatives of the singularity 
potential in combination with Eqs. (3.44a), (12a), (12b), (12c), and (23c). This solution for 
the upwash velocity showed good agreement with a classical solution given by Eq. (2.30) 
in Ref. [2]. 

Summary 

Singularities like point sources, point doublets, vortices, and line doublets may be used 
to compute wind tunnel wall interference effects. These singularities represent blockage 
and lifting effects of a wind tunnel model. Their location and strength is closely related to 
the geometry of the test article. Therefore, it is necessary to apply the Prandtl-Glauert 
transformation to the singularity strength as well, if wall interference effects are computed 
in a subsonic flow field. The table below summarizes the particular form of the Prandtl- 
Glauert transformation for each singularity type that is discussed in this section. 

Table 7 : Singularity Strength Transformation {o' = a ■ F((3) ) 


Singularity Type 

Singularity Strength (a) 

FW 

Point Source (2D) 

Uqq • D 

8 

Point Doublet (2D) 

Uoo ■ 2 • 7T • r\ 

8 

Point Source (3D) 

Uoo ■ 7T • rlc 

8 2 

Point Doublet (3D) 

Uoo ■ 2 • 7T • rl 

8 2 

Vortex (2D) 

L / {poo ‘ u oo * As) — F 

8 2 

Line Doublet (3D) 

L / {poo * ^ oo) = F * As 

8* 


Important Remark : Previously, the derivation of the singularity strength transforma- 
tion for a three-dimensioned line doublet was performed erroneously by the author (Eq. (22) 
in Ref. [1] is not valid for a line doublet). The correct singularity strength transformation 
for a line doublet is now given in Eq. (3.44a). The author added a corresponding errata 
page to Ref. [1]. 
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APPENDIX 4 


METHOD OF IMAGES - POINT DOUBLET 


The Method of Images may be applied to compute the wind tunnel and wall inter- 
ference flow field caused by a point doublet that is located inside of a rectangular wind 
tunnel with “closed wall” or “open jet” boundary conditions. Figure 16a shows the image 
system for the “closed wall” boundary condition. All unit point doublet strength values 
have the same sign in this case. Figure 16b shows the image system for the “open jet” 
boundary condition. Point doublet strength sign changes are indicated in Fig. 16b. 

It is assumed that a point doublet points in the negative x-direction and is located 
at x 0 ,y 0 ,z 0 inside of a rectangular tunnel of width hi and height h 2 (see Fig. 16a). Then, 
we get for the velocity components u = d<j>j dx, v = d<f>/dy , w = d<p / dz of the wind tunnel 
flow field at a flow field point x,y,z : 


where 


OO OO 


II 

H 

E 

m = — oo 

U (m, n) - F(m, n ) 

n— — oo 

(4.1a) 

II 

N 

oo 

E 

m= — oo 

oo 

^ V(m,n) • F(m,n) 

n = — oo 

(4.16) 

w(x,y,z) = 

oo 

E 

oo 

y, W(m,n) ■ F(m,n) 

(4.1c) 


m — — oo n= — oo 


_f_ {Y{m)-y o) 2 + (Z(n)-x 0 ) 2 - 2-(x-x 0 ) 2 x 

4 7 T ~f ‘ TTfT (4.2a) 

(x - x 0 ) 2 + ( Y(m ) - y 0 ) 2 + (Z(n) - z 0 ) 2 


V(m, n) 
W(m, n) 




and 


y(m) 

Z(n) 


y + m • h\ + [ 1 — (— l) m ] ■ Vo 


(4.3 a) 
(4-3 b) 


z + 


n • /12 + [1 — (—1)" ] ’ 


The velocity components u w — d<i> w /dx, v w = d<f> w /dy , w w = d(j> w ldz of the wall 
interference flow field at a flow field point x,y,z may be obtained by only considering 
summation terms in Eqs. (4.1a), (4.1b), (4.1c) that are caused by image point doublets. 
Then, we get : 

OO OO 

E E U(m,n ) • F(m, n ) 


m= — oo n=— oo 
(m,n)^(0 t 0) 
oo oo 


,(x,y,z) = 52 V(m,n)-F(m,n) 


(4.4a) 


(4.46) 


(4.4c) 


to— — oo n= — oo 

(m,n)7*(0,0) 

oo oo 

w u ,(x,y,z) = 52 W( m .rc)--F(m,n) 

rn=— oo n= — oo 
(m,n)?£(0,0) 

Perturbation velocity components £7(m,n), V(m,n), W(m,n), are given by Eqs. (4.2a), 
(4.2b), and (4.2c) . Factor F(m, n) is given by Eq. (4.2d) . 

Numerical studies have shown that summation terms in Eqs. (4.1a), (4.1b), (4.1c), 


(4.4a), (4.4b), an d (4.4c) should be computed using “DOUBLE PRECISION” arithmetic 
(“*r8” compiler flag). A summation index range of —200 <m< 200 and —200 < n < 200 
is recommended. 
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APPENDIX 5 


METHOD OF IMAGES - LINE DOUBLET 


The Method of Images may be applied to compute the wind tunnel and wall interfer- 
ence flow field caused by a semi-infinite line doublet that is located inside of a rectangular 
wind tunnel with “closed wall” boundary condition. The Method of Images adds pertur- 
bation velocity components of each line doublet of the image system to obtain a wind 
tunnel and wall interference flow field solution (see Fig. 17a). Unfortunately, perturbation 
velocity components of a line doublet are a function of the line doublet orientation angle 
0. This angle is defined as 0 [ rad] if the lift force caused by the line doublet points in 
the positive z-axis direction of the tunnel coordinate system x,y, z (see Fig. 17a). A co- 
ordinate system transformation may be used to simplify the application of the Method of 
Images to a line doublet (for more detail see App. 2) . This transformation is essentially 
a coordinate system rotation that uses line doublet orientation angle 0. 

It is assumed that a semi-infinite line doublet with starting point x 0 ,y 0 ,z 0 and orien- 
tation angle 0 is located inside of a rectangular wind tunnel of width hi and height hi 
(see Fig. 17a). Then, based on the Method of Images and Eqs. (2.10a), (2.10b), (2.10c) in 
App. 2, we get for velocity components u — d<f>/dx, v = d<f>/dy, w = d<f>/dz of the wind 
tunnel flow field at a flow field point x,y,z : 

OO OO 

-<*.»■*)- E E U(m,n) (5.1a) 

m= — oo n= — oo 


OO oo 

v(x,y,z) = ^2 

m= — o o 7i = — oo 
oo oo 

™(x,y, z) = ^ ^ [ -F(m,n)'Sm[r(m,n)] + W(m,n) • cos[T{m,n)\ 

m = — oo n = — oo 


V(m,n) ■ cos[T(m,n)] + W{m, n) ■ sm[r(m, n)] 


(5.16) 


(5.1c) 


where velocity components U(m,n), V(m,n), and W(m,n ) are computed in the line 
doublet fixed coordinate system of each line doublet of the image system. Using 

Eqs. (2.4), (2.5a), (2.5b), (2.5c), (2.6a), (2.6b), (2.6c), from App. 2, line doublet velocity 
components may be written as : 

° C - Co 


Z7(m,n) = 


4 7 r 


(C ~ Co) 2 +(r) - rj 0 ) 2 + (C - Co) 2 


1 3/2 


(5.2) 
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V(m,n) = 


where 


— a 


iv - rio) ■ (C - (o) 

4 7T (r] - TJo) 2 + (C - Co) 2 


2 B 


(V~Vo) 2 + (C - Co) 2 


+ c 


1 + 


C"Co 


c = 


« - {.)* + (I? — l.) 2 + (C - c .) 2 ] 
( - u 


1/2 


(C - C o ) 2 + (v - Vo) 2 + (C - Co) 2 


3/2 


W(m,n) = 


where 


4 7T (*7 — Vo) 2 + (C - Co) 2 


D = 1 + 


(q-qo) 2 -(C-Co) 2 
0? - »7o) 2 + (C - Co) 2 

C - Co 


•D + E 


E = 


(C - Co) 2 + (v~ Vo) 2 + (C - Co) 2 ] 
(-1) • (C - Co) • (C - Co) 2 


1/2 


-, 3/2 


(5.3a) 


(5.36) 


(5.3c) 


(5.4a) 


(5.46) 


(5.4c) 


(( - to ) 2 + (V ~ Vo ) 2 + (C ~ to ) 2 

The coordinates of the flow field point x , y, 2 in the tunnel coordinate system are 
transformed to the line doublet fixed coordinate system £, Vi C using the following equations 
(see also Eqs. (2.1a), (2.1b), (2.1c)) : 


f = x (5.5a) 

77 = 77 ( 771 , n) = Y(m) • cos[r(m,7i)] - Z(n) ■ sin[r(m, n)] (5.56) 

( = £( 771 , 77 ) = K(m)-sin[r(m,n)] + Z(n) • cos[r(m, 71 )] (5.5c) 

The Y and Z coordinate of a flow field point have to be computed using the following 
equations : 

Y(m) = y + m • h\ + [ 1 — (— l) m ] • y 0 (5.6a) 

Z(n) = z + n-h 2 + (1 - (-1)" ]-Zo (5.66) 


Similarly, coordinates of the line doublet starting point r o ,yo,z 0 are transformed to 
the line doublet fixed coordinate system tiVit using the following equations (see also 
Eqs. (2.2a), (2.2b), (2.2c)) : 

Co = (5-7a) 
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J7 0 (m,n) = j/o • cos[r(m, n)] — z 0 • sin[r(m, n)] (5.76) 

Co(m,n) = y 0 • sin[r(m, n)j + z Q ■ cos[r(m, n)] (5.7 c) 

The line doublet orientation angle of each line doublet of the image system is a func- 
tion of line doublet orientation angle 0 and image system indices (m, n). By inspection 
(Fig. 17a), we see that the line doublet orientation angle of each image system element for 
the “closed wall” boundary condition is given as : 

r(m,n) = (-l) m • [ (-l) n • 0 + n • x ] ; closed wall (5.8) 


where 0 = t( 0,0) is the angle between the z— axis of the tunnel coordinate system x,y,z 
and the ((-axis of the line doublet fixed coordinate system £, rj, ( . 

The velocity components u w = d<j) w jdx , v w = d<f> w /dy, w w = d<t> w /dz of the wall 
interference flow field at a flow field point x,y,z may be obtained by only considering 
velocity components in Eqs. (5.1a), (5.1b), (5.1c) that are caused by image line doublets. 
Then, we get : 

oo oo 


u w (x,y,z ) = ^ U(m,n) 



m= — oc n= — oo 

(5.9a) 


(m,n)^(0,0) 


OO 

oo 


<2 

6 

II 

M 

y [ V(m,n) • cos[T(m, n)j + W(m,n) ■ sm[r(m,n)] ] 


J 

m — — oo n=— oc 

(5.96) 

(m,n 

)^(0,0) 


oo 

OO 


w w (x,y,z) = ^ 

y — V(m,n) ■ sin[r(m,n)} + W(m,n) ■ cos[T(m,n)\ 

L 


771= — OO 

n= — oo 

(5.9c) 


(m,n)?£(0,0) 


Velocity components £/(m,n), V(m,n), TV(m,n), and orientation angle r(m,n) axe 
computed using Eqs. (5.2), (5.3a), (5.3b), (5.3c), (5.4a), (5.4b), (5.4c), and (5.8) . 

Numerical studies have shown that summation terms in Eqs. (5.1a), (5.1b), (5.1c), 
(5.9a), (5.9b), and (5.9c) should be computed using “DOUBLE PRECISION” arithmetic 
(“-r8” compiler flag). A summation index range of -200 < m < 200 and -200 < n < 200 
is recommended. 

The Method of Images may also be used to compute the wall interference flow field of 
a line doublet located in a rectangular “open jet” if the line doublet orientation angle 0 
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is equal to 0, tt/2, tt, or 3 tt/ 2. Figure 17b shows the image system of a rectangular “open 
jet” for 0 = 0 . The orientation angle of each image system element of the “open jet” 
becomes : 


r(m, n) 


J 0 + m ■ 7 r for 0 = 0 or 7r 

^ 0 + n • 7 r /or 0 = 7 t/ 2 or 37 t/2 


open jet (5.10) 


Thus, Eq. (5.8) has to be replaced by Eq. (5.10) if the interference flow field of a line doublet 
in an open jet for any of the four possible values of 0 is computed by using Eqs. (5.9a), 
(5.9b), and (5.9c) . 
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PART 2 

SOLUTION OF LINEAR SYSTEM 

(LU - DECOMPOSITION ALGORITHM) 

I I 

PART 3 

CALCULATION OF WALL INTERFERENCE FLOW FIELD 

1. SUMMATION OF WALL PANEL CONTRIBUTIONS AT FLOW FIELD POINTS 
2. REVERSAL OF PRANDTL - GLAUERT TRANSFORMATION 



Fig. 2 Basic structure of panel method code ANTARES. 
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Fig. 4 Panel representation of square wind tunnel. 







CLOSED WALL BOUNDARY CONDITION ; FULLSPAN CONFIGURATION 
UNIT POINT DOUBLET AT X = 0 [FT] , Y = -1 [FT] , Z = -2 [FT 
WALL INTERFERENCE ON LINE 1 (Y = 0 [FT] , Z = 0 [FT]) 



0.0006 





CLOSED WALL BOUNDARY CONDITION ; FULLSPAN CONFIGURATION ; MACH 
UNIT POINT DOUBLET AT X = 0 [FT] , Y = -1 [FT] , Z = -2 [FT] 
WALL SIGNATURE ON ROW 1 (Y = 5 [FT] . Z = 2 [FT]) AND ROW 2 (Y = -5 [FT] 



Fig. 5e Velocity component u/u^ (Moo = 0.80) due to unit point doublet in 
square tunnel; Method of Images (M.I.) versus panel method code (P.M.). 





CLOSED WALL BOUNDARY CONDITION ; FULLSPAN CONFIGURATION 
UNIT LINE DOUBLET AT X = 0 [FT] , Y = -1 [FT] . Z = -2 [FT] ; ORIENT. ANGLE 
WALL INTERFERENCE ON LINE 1 (Y = 0 [FT] , Z = 0 [FT]) 



Fig. 6a Wall interference velocity component u^/itoo due to unit line doublet in 
square tunnel; Method of Images (M.I.) versus panel method code (P.M.). 



CLOSED WALL BOUNDARY CONDITION ; FULLSPAN CONFIGURATION 
UNIT LINE DOUBLET AT X = 0 [FT] , Y = - 1 [FT] , Z = -2 [FT] ; ORIENT. ANGLE 
WALL INTERFERENCE ON LINE 1 (Y = 0 [FT] , Z = 0 [FT]) 


O 



Fig. 6b Wall interference velocity component v w /u^ due to unit line doublet in 
square tunnel; Method of Images (M.I.) versus panel method code (P.M.). 


0.0006 



CLOSED WALL BOUNDARY CONDITION ; FULLSPAN CONFIGURATION 
UNIT LINE DOUBLET AT X = 0 [FT] , Y = -1 [FT] , Z = -2 [FT] ; ORIENT. ANGLE 
WALL INTERFERENCE ON LINE 1 (Y = 0 [FT] , Z = 0 [FT]) 



Fig. 6c Wall interference velocity component i0ti>/woo due to unit line doublet in 
square t unn el; Method of Images (M.I.) versus panel method code (P.M.). 




CLOSED WALL BOUNDARY CONDITION ; FULLSPAN CONFIGURATION ; MACH = O.C 
UNIT LINE DOUBLET AT X = 0 [FT] , Y = -1 [FT] , Z = -2 [FT] ; ORIENT. ANGLE 
WALL SIGNATURE ON ROW 1 (Y = 5 [FT] , Z = 2 [FT]) AND ROW 2 (Y = -5 [FT] , Z 



Fig. 6d Velocity component u/u 0 c (Moo — 0.00) due to unit line doublet in 
square tunnel; Method of Images (M.I.) versus panel method code (P.M.). 



CLOSED WALL BOUNDARY CONDITION ; FULLSPAN CONFIGURATION ; MACH = 0.8 
UNIT LINE DOUBLET AT X = 0 [FT] , Y = -1 [FT] , Z = -2 [FT] ; ORIENT. ANGLE 
WALL SIGNATURE ON ROW 1 (Y = 5 [FT] , Z = 2 [FT]) AND ROW 2 (Y = -5 [FT] , Z = 

□ M.I.(ROW 1) O M.I.(ROW 2) * P.M.(ROW 1) A P.M.(ROW 2) 



Fig. 6e Velocity component u/u TO (Moo = 0.80) due to unit line doublet in 
square tunnel; Method of Images (M.I.) versus panel method code (P.M.). 
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ROW 2 I ^ ROW 1 
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Fig. 7 Geometry of rectangular wind tunnel and location of singulariti 
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Fig. 8 Panel representation of rectangular wind tunnel. 
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Fig- Wall interference velocity component u w ju 0 0 due to two unit point doublets 
rectangular tunnel; Method of Images (M.I.) versus panel method code (P.M.)- 
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Fig. 9b Wall interference velocity component v w /u 00 due to two unit point doublets m 
rectangular tunnel; Method of Images (M.I.) versus panel method code (P.M.). 






CLOSED WALL BOUNDARY CONDITION ; SEMISPAN CONFIGURATION ; MACH = 0.00 
UNIT POINT DOUBLET AT X = 0 [FT] , Y = -1 [FT] , Z = -2 [FT] AND X = 0 [FT] , Y = -1 [FI 
WALL SIGNATURE ON ROW 1 (Y = 5 [FT] , Z = 2 [FT]) AND ROW 2 (Y = -5 [FT] , Z = 



Fig. 9d Velocity component u/uqc (Moo = 0.00) due to two unit point doublets in 
rectangular tunnel; Method of Images (M.I.) versus panel method code (P.M.). 



CLOSED WALL BOUNDARY CONDITION ; SEMISPAN CONFIGURATION ; MACH = 0.80 
UNIT POINT DOUBLET AT X = 0 [FT] , Y = - 1 [FT] , Z = -2 [FT] AND X = 0 [FT] , Y = -1 [FT] 
WALL SIGNATURE ON ROW 1 (Y = 5 [FT] , Z = 2 [FT]) AND ROW 2 (Y = -5 [FT] , Z = ; 


INI Li_ 



Fig. 9e Velocity component u/u Q 0 (M 0 c = 0.80) due to two unit point doublets in 
rectangular tunnel; Method of Images (M.I.) versus panel method code (P.M.). 



CLOSED WALL BOUNDARY CONDITION ; SEMISPAN CONFIGURATION 



Fig. 10a Wall interference velocity component u w juca due to two unit line doublets in 
rectangular tunnel; Method of Images (M.I.) versus panel method code (P.M.). 


0.004 






'SED WALL BOUNDARY CONDITION ; SEMISPAN CONFIGURATION 
0 [FT] , Y = -1 [FT] , Z = -2 [FT] AND X = 0 [FT] , Y = 


Fig. 10c Wall interference velocity component w w /uoo due to two unit line doublets in 
rectangular tunnel; Method of Images (M.I.) versus panel method code (P.M.). 











CLOSED WALL BOUNDARY CONDITION ; SEMISPAN CONFIGURATION ; MACH = 0.80 
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Fig. 11a Wall interference velocity component u w / due to unit point doublet in 
open jet; Method of Images (M.I.) versus panel method code (P.M.). 
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OPEN JET BOUNDARY CONDITION ; OPEN JET DIMENSION^ : 1 
UNIT POINT DOUBLET AT X = 0 [FT] . Y = -1 [FT] , Z 
WALL INTERFERENCE ON OPEN JET CENTERLINE (Y = 0 [FT 



pig. nb Wall interference velocity component v w /uoo due to unit point doublet in 
open jet; Method of Images (M.I.) versus panel method code (P.M.)- 
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OPEN JET BOUNDARY CONDITION ; OPEN JET DIMENSIONS : 10 [FT 
UNIT LINE DOUBLET AT X = 0 [FT] , Y = -1 [FT] , Z = -2 [FT] ; OR 
WALL INTERFERENCE ON OPEN JET CENTERLINE (Y = 0 [FT] , Z 

□ M.I.(MACH=0.00) O M.I.(MACH=0.90) * P.M.(MACH=0.00) A I 



Fig. 12a Wall interference velocity component u w /uoo due to unit line doublet in 
open jet; Method of Images (M.I.) versus panel method code (P.M.). 


0.0010 





°°n/ m a 

Fig. 12b Wall interference velocity component u^/uoo due to unit line doublet in 
open jet; Method of Images (M.I.) versus panel method code (P.M.). 
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Fig. 12c Wall interference velocity component w w /uoo due to unit line doublet in 
open jet; Method of Images (M.I.) versus panel method code (P.M.). 
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Fig. 13 Blockage factor (u w /u 0 Q ) on tunnel centerline due to 2D unit point doublet 
in 2D wind tunnel; NACA TN 3176 (Exact Solution) versus ANTARES (Panel Method). 
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Fig. 14 Upwash velocity (u> w /uoo) on tunnel centerline due ^^D^ut vortex 

in 2D wind tunnel; AEDC-TR-69-47 (Exact Solution) versus ANTARES (Panel Method). 
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Fig. 15 b Perturbation velocity component transformation. 
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METHOD OF IMAGES - POINT DOUBLET 
(OPEN JET BOUNDARY CONDITION) 
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Fig. 16 b Image system of point doublet for open jet wall boundary condition. 
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METHOD OF IMAGES - LINE DOUBLET 
(CLOSED WALL BOUNDARY CONDITION) 
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METHOD OF IMAGES - UNE DOUBLET 
(OPEN JET BOUNDARY CONDITION) 
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Fig. 17 b Image system of line doublet for open jet wall boundary condition. 
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